Vison gap in the Rokhsar-Kivelson dimer model on the triangular lattice 
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With the classical Monte Carlo method, I find the energy gap in the Rokhsar-Kivelson dimer 
model on the triangular lattice. I identity the lowest excitations as visons, and compute their energy 
as a function of the momentum. 



The Rokhsar-Kivelson (RK) dimer model is an in- 
teresting example of a two-dimensional quantum liquid 
[1-3] . When the kinetic and potential coupling constants 
of the RK dimer model are equal (the so-called "RK 
point"), the ground state of the model is known exactly, 
and the ground-state properties on planar lattices may be 
analytically studied with the Pfaffian method [4]. Much 
attention was paid to the study of the RK model on the 
triangular lattice, in which case the RK point is known 
to be in the liquid phase with topological order and with 
exponentially decaying correlations [2,5,6] (on the square 
lattice, in contrast, the RK point has power-law correla- 
tions and a gapless spectrum, due to the bipartite nature 
of the lattice [1,3,7-9]). Unfortunately, our knowledge of 
the excitation spectrum at the RK point on the triangu- 
lar lattice is limited to numerical methods [2,10]. Study- 
ing excitations in the RK model is interesting in view 
of the prediction that the elementary excitations are the 
so-called "visons" which are non-local objects in terms of 
dimers. Visons were proposed as elementary excitations 
in quantum liquids with gauge symmetry (possibly 
also in spin- 1/2 systems) [11,12]. While applicability of 
this proposal to realistic spin- 1/2 systems is debatable, 
there exist model systems, where the existence of vison 
excitations is explicitly shown and the whole spectrum of 
visons may be exactly found [13,14]. The usual price to 
pay for the exact solvability is that in such models the 
vison excitations are strictly local (i.e., their energy is 
independent of the momentum) and non-interacting. In 
contrast, the RK model presents a very non-trivial case, 
where the visons have their dispersion and interact in 
many-vison states. 

In this work, I use the proposal of Henley to extract the 
energy gap of the excitations from the time correlation 
of the classical Monte Carlo algorithm used to compute 
the ground-state observables [9] (the classical calculation 
of the excitation spectrum is actually possible due to the 
supersymmetry of the RK point [15]; the relation be- 
tween supersymmetric quantum systems and stochastic 
classical dynamics was discussed in various contexts [16]). 
By choosing appropriate correlation functions, we can re- 
solve the momentum-dependence of the gap, and also dis- 
tinguish excitations carrying odd number of visons from 
those carrying even number of visons. We find that the 
low-energy excitations are indeed vison-like (carry odd 
number of visons) and plot their dispersion in the Bril- 



louin zone. 

The Rokhsar-Kivelson dimer model may be defined on 
any graph: the dimer coverings of the graph define an or- 
thonormal basis of the Hilbert space; and the quantum 
Hamiltonian is 
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RK 
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The sum is taken over all length-four loops of the 
graph. The two coupling constants t and v determine the 
strength of kinetic and potential terms, respectively. At 
the "RK point" t = v, the Hamiltonian can be shown to 
be non- negative (assuming t > 0), and its ground state 
may be constructed as the sum of all dimer configura- 
tions taken with equal amplitudes (more precisely, we 
can restrict the sum to the dimer configurations from 
any connected component of the "phase space" , i.e. to 
configurations which can be obtained from each other 
by the kinetic term in the Hamiltonian (1)) [1]. This 
ground state has energy zero, which is a manifestation 
of the supersymmetry of the RK point [15]. We further 
specify to the case of the underlying graph being the two- 
dimensional triangular lattice (the sum in the Hamilto- 
nian (1) is taken over all rhombi, so the total number 
of terms in the sum is three times the number of lattice 
sites), and set the energy units t = v = 1. Our main ob- 
jective in this paper is finding the low-lying excitations 
of this model. From earlier numerical studies, it has been 
suggested that this model has a gapped spectrum [2,10], 
which seems to be in agreement with the exponential 
decay of ground-state correlation functions known from 
analytic studies [2,5,6]. 

From the general discussions of the dimer liquids in 
two dimensions, the Z%- vortex operator (the so called "vi- 
son" ) is known to play an important role among physical 
observables. The vison (more precisely, the "two- vison") 
operator Vr is defined for any contour T intersecting links 
of the lattice and terminating either at the lattice bound- 
ary or inside a plaquet in the bulk of the lattice (Fig. la). 
The operator Vr is defined as the parity of the number 
of dimers intersecting T: 



V r = (-1 



,no. of dimers intersecting T 



(2) 



If one commutes such an operator with the Hamiltonian 
(1), the only non- vanishing contribution comes from the 
rhombi containing the end points of the contour. In par- 
ticular, if the contour forms a closed loop or terminates 
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at a lattice boundary, the corresponding operator exactly 
commutes with the Hamiltonian and gives rise to differ- 
ent topological sectors of the Hilbert space [1,11,17]. 




FIG. 1. (a) Definition of the vison operator Vr- The con- 
tour T (dashed line) connects two triangular plaquets of the 
lattice (shaded). The value of Vr is the parity of the intersec- 
tion of dimers with V. For the contour F and for the dimer 
configuration shown in the figure, Vr = — 1 (three intersec- 
tions), (b) The reference dimer configuration for the sign 
fixing of a single- vison operator. An alternative (to that de- 
scribed in the main text) formulation of the sign-fixing rule: 
the contours T must be drawn in such a way that they do not 
intersect dimers from the reference configuration. 

It is natural to suggest that terminating the contour 
in the bulk of the lattice produces an excited state close 
to an cigenstate. Of course, to form a true excited eigen- 
state would require some "dressing" of the vison opera- 
tor with local corrections around the contour end point. 
However, the topological structure of the excitation will 
be preserved vison-likc. To clarify the above reasoning, 
we first describe some properties of the vison operators 
Vr. First, it is easy to check that the operator Vr, up 
to a sign, depends only on the end points of the contour; 
the dependence on the contour itself reduces to a control- 
lable change of the sign: changing from the contour T to 
another contour T' with the same end points changes the 
sign of the vison operator as Vr = (— 1) Vr', where S is 
the number of lattice points between the contours T and 
r'. Second, concatenation of the contours corresponds 
to the multiplication of the corresponding vison opera- 
tors. Therefore we may represent the operator Vr as the 
product of two vison operators at end points: Vr = V1V2, 
where each of the "single-vison" operators Vi and V2 de- 
pends on one of the two end points of T. Constructed 
this way, the point vison operators V, obey Zi algebra 
(V^ 2 = 1) and are defined on the frustrated dual lattice: 
their index i refers to a plaquet of the original lattice, 
and they change sign on going around one lattice site of 
the original lattice. For the triangular lattice, we should 
think of visons as living on the hexagonal lattice with the 
magnetic flux of half quantum per hexagon. 

To establish a sign convention for visons, we need to 
fix a Z2 gauge on the dual lattice. This can be most 
easily done by taking a certain (arbitrary, but fixed once 
forever) dimer configuration as a reference one. Then 
in the definition (2), we multiply the right-hand side by 
the same expression Vr computed in the reference dimer 
configuration. With the new definition, Vr becomes a 



single-valued function of the end points of T (indepen- 
dent of the choice of the contour). In our calculation, 
we take the reference dimer configuration as shown in 
Fig. lb (note an additional doubling of the unit cell). 

An important property of the vison operator is that it 
is a non-local operator in terms of dimers. A single-vison 
operator V. in an infinite system involves the contour T 
continued to infinity and hence corresponds to a change 
in the boundary conditions on the wave function at in- 
finity. Namely, a circular permutation of dimers along a 
big closed contour encircling the "excitation region" re- 
verses sign of the excitation with odd number of visons, 
but keeps the wave function unchanged for the excitation 
with even number of visons (Fig. 2). According to this 
criterion, we may classify any "localized" wave packet as 
either vison-like or non-vison-like. Non-vison-like excita- 
tions are expressed as local operators in terms of dimers. 
A vison-like excitation may be described as a non-vison- 
like (local) operator multiplied by a vison operator. Thus 
we naturally have two classes of excitations with a Z2 
grading: combining two vison-like excitations we obtain 
a non- vison like excitation, while combining a non-vison- 
like excitation with a vison-like excitation gives a vison- 
like excitation. Of course, with this construction it seems 
natural that vison-like excitations should be considered 
as "elementary" excitations, while non-vison-like excita- 
tions may be constructed as composite excitations from 
vison-like elementary excitations. In reality, however, it 
may happen that vison-like excitations are pushed high 
in energy, so that the low-lying physical excitations are 
all non-vison-like. In this paper I demonstrate that it is 
not the case at the RK point on the triangular lattice. 
We shall see that both the vison-like and non-vison-like 
sectors are gapped and that the gap in the vison-like sec- 
tor is smaller than in the non-vison-like sector. So vison 
excitations are indeed the lowest-energy excitations in 
the model. 




vison cut 



FIG. 2. A localized wave packet of vison-like excita- 
tions. Vison-like excitations may be distinguished from 
non-vison-like ones by a circular permutation of dimers along 
a big contour (shown). If we denote such a permutation by 
P, vison-like excitations obey the equation P 2 = — P, while 
non-vison-like excitations have P 2 = P. 

The exponentially decaying ground-state correlation 
functions [2,5] suggest a gap in the excitation spectrum, 
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and indeed both quantum Monte Carlo studies [2] and 
exact diagonalization on small systems [10] suggest the 
presence of the gap (from exact diagonalization, the value 
of the gap was estimated as 0.1). However, as pointed 
out by Henley [9] , at the RK point the gap may be more 
easily extracted from a classical Monte Carlo simulation 
similar to that used for calculating the ground-state ex- 
pectation values (see, e.g., Rcfs. [18,6]). Namely, con- 
sider the following random walk defined on the space of 
all dimer coverings of the lattice. A step of the random 
walk is defined as picking at random any rhombus and, 
if it contains two parallel dimers, flipping this pair of 
dimers into the other two sides of the rhombus (as indi- 
cated by the kinetic term of the Hamiltonian (1)). If the 
chosen rhombus is non-fiippable, the dimer configuration 
remains unchanged at this step of the random walk. 

As shown by Henley [9], such a random walk sim- 
ulates the quantum-mechanical evolution in imaginary 
time. Accordingly, the exponents governing the decay 
of the dynamic correlation functions with time are given 
precisely by the excitation energies of the quantum sys- 
tem. For our discrete random walk, this procedure of 
determining the excitation energies produces systematic 
errors arising from the discretization of time steps. The 
discreteness of time steps may, in principle, be prop- 
erly compensated; however, we simply neglect the cor- 
responding systematic errors. One can show that dis- 
cretization of time leads to relative corrections to the gap 
magnitude of approximately one over the total number of 
rhombi in the system. In our Monte Carlo calculation we 
take a sufficiently large system of 20x20 sites (thus con- 
taining 1200 rhombi), and those corrections are smaller 
than the statistical errors for the lengths of random walks 
used in our calculations. Therefore we disregard the 
time-discretization errors and extract the gaps directly 
from the correlations of the discrete random walk. 

The next useful observation is that, by taking appro- 
priate correlation functions, we can probe the gap at a 
given wave vector. And, moreover, we can distinguish 
between vison-like and non-vison-like excitation sectors. 
For computing the gap in the non-vison-like sector, we 
should consider correlations of non-vison-like observables 
(e.g., the dimer density). For the gap in the vison sector 
we take vison-like observables (e.g., the point vison Vi 
defined above). 

We first consider the excitations in the vison-like sec- 
tor. In order to compute the gap, we take the correlation 
function 

F(v ij ,t-t') = (V i (t)V j (t')), (3) 

where Vi(t) is the point vison on the triangle i at the mo- 
ment t of the random-walk procedure as defined above, 
Yij is the vector connecting the plaquets i and j. This 
correlation function is properly defined, in spite of hav- 
ing only one vison operator at each of the time mo- 
ments t and t' . To demonstrate the consistency of 



the definition (3), it is convenient to insert the square 
of the vison operator [V^(t)] 2 = 1. Then we rewrite 
Vi(t)Vj(t') = {VityVjityiVji^Vji?)). The first product 
of the two vison operators is defined in (2) with the con- 
tour r connecting the points i and j. The second product 
involves two vison operators at one space point, but at 
different time moments. Obviously this product is also 
well defined for our random-walk process: namely, every 
dimer flip at any time between t and t' on a rhombus con- 
taining the triangle j changes the sign of Vj(t)Vj(t'). In 
other words, Vj(t)Vj(t') counts the parity of the number 
of such flips. 




FIG. 3. This figure shows the Brillouin zones for 
non-vison-like and vison-like excitations (superimposed, on 
the same scale). Non-vison-like excitations: the Brillouin 
zone is the big hexagon (solid line; the side length of the 
hexagon is 4n/3 in the units of the inverse lattice constant of 
the original triangular lattice) , with the minimum-gap points 
marked with big circles and labeled a. Vison-like excita- 
tions: the Brillouin zone is the dashed rectangle (for the 
choice of gauge specified in Fig. lb). The small solid circles, 
squares, and triangles (labeled with the letters A, B, and 
C, respectively) are the high-symmetry points. For exam- 
ple, points A are centers of six-fold symmetry. The vison gap 
is found to reach its minimum at points B, and to increase 
strongly towards points A. 

With the gauge choice for visons as discussed above 
(Fig. lb), the unit cell of the lattice is doubled, and 
contains four triangular plaquets. To characterize vi- 
son excitations with wave vectors, we perform a Fourier 
transform of the correlation function (3) and arrive at 
the 4x4 matrix F(k, t), where k is the wave vector from 
the reduced Brillouin zone. The full Brillouin zone and 
the reduced Brillouin zone are sketched in Fig. 3. Note 
that with our choice of the vison gauge, the correlation 
function F(k, t) has certain symmetries in the fc-space re- 
flecting the symmetries of the original triangular lattice. 
Those symmetries contain the six- fold rotation/reflection 
symmetry (D 6 symmetry group) around a certain point 
in the fc-space, together with translations by two basis 
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vectors defining a triangular lattice. The high-symmetry 
points in the /c-space are marked in Fig. 3 by letters A, 
B, and C. 
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FIG. 4. Top: the energy of the vison-like excitations along 
the section A-B-C-A of the Brillouin zone shown in Fig. 3. 
The error bars are much smaller than the symbol size. The 
arrows on the energy axis indicate positions of As (bottom of 
the vison band), Ao (lowest-energy non- vison-like excitation), 
and Ab + Ao (an estimate for the bottom of the three- vison 
continuum). Bottom: three-dimensional plot of the energy 
of the vison-like excitations as a function of the wave vector. 
The center and the corners of the hexagon correspond to the 
points A, as labeled in Fig. 3. 

Now we determine the energy of the vison excitations 
A(k) from the exponential decay of the diagonal ele- 
ments of F(k,i) [19]. The sections of the energy disper- 
sion along the "crystallographic axes" of the Brillouin 
zone and the 3-D plot of A(k) are presented in Fig. 4. 
The points B of the fc-space give the minimal energy gap 
Ab = 0.089(1) [which agrees with Ref. [10] claiming the 
gap value 0.1]. The points C are the saddle points of 
the energy dispersion with Ac = 0.114(1). The points A 
are the centers of high-energy regions. In those regions, 
a naive fitting with an exponential suggests A^ > 0.3 
and does not reproduce well the i-dependence of F(k, t). 
This may indicate that, at those wave vectors, the lowest 
excitations are not elementary visons, but rather combi- 
nations of three visons from the points B. Thus the t- 
dependence of F(k, t) reflects not an isolated excitation, 
but the bottom of a multi-particle continuum. The dif- 



ference in the time dependence between low-energy and 
high-energy regions in the k space is illustrated in Fig. 5 
where we show typical i-dependences of F(k, t) for the 
two /c-points: one near point B in the fc-space, and the 
other one near point A. 
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FIG. 5. Typical vison and dimer correlations in the clas- 
sical Monte Carlo random walk. Solid an dashed line show 
the vison correlations for wave vectors near points B and A, 
respectively (see Fig. 3). Note that while near point B the 
correlation function has a well-marked exponential behavior 
at long times, the correlations near point A do not allow a 
precise fitting with a single exponent. The dotted line shows 
the dimer-dimer correlation function at the point a (the mini- 
mum-gap point). A good exponential decay of the correlation 
function possibly indicates a two- vison bound state. 

Next, we repeat the same procedure for non-vison- 
like excitations by taking the correlations of the dimer- 
density operator instead of the vison operator Vi{t) in (3). 
We then find that the lowest gap in the non-vison sector 
is Ao = 0.144(1) and is reached at points labeled a in 
Fig. 3. The t-dependence of the dimer-dimer correlation 
function at point a is shown in Fig. 5. A good fit with 
an exponential dependence, together with the inequality 
Ao < 2 As, suggests that the lowest non- vison-like exci- 
tation is not just a superposition of two non-interacting 
vison excitations, but a bound state of such a pair. On 
the other hand, its energy is considerably higher than 
that of an elementary vison excitation (Ab), which con- 
firms the claim that the lowest excitation is vison-like. 
[This result also suggests to estimate the bottom of the 
continuum in the vison-like excitations as A + Ab in- 
stead of 3A^, see Fig. 4a.] 

Finally, with this method of calculating the excitation 
gap, we can verify the claim of Ref. [10] about the absence 
of low- lying edge states in the case of lattices with bound- 
aries. We examine the excitation spectrum of the 10x10 
and 20x20 cylinders with straight boundaries along the 
lattice directions. From the result on non-vison-like ex- 
citations in the bulk, we expect that visons are attracted 
to each other, and therefore, visons should also get at- 
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tracted to boundaries (since the vison cut may be ter- 
minated at the boundary at no energy cost). Note also 
that near the boundary there is no distinction between 
vison-like and non- vison- like excitations. Thus for deter- 
mining the energy of the edge states we may take the 
dimer-dimer correlation function at the very boundary. 
Our classical Monte Carlo simulation gives the boundary 
gap Acdgo = 0.072(1) reached at the wave vector n along 
the boundary. In line with our expectations, the gap at 
the boundary is indeed somewhat reduced compared to 
the bulk vison gap, but remains finite. 

The author thanks M. Feigelman for helpful discus- 
sions and comments. 
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